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Abstract 



Let (M, g) be a compact Riemmanian manifold with non-empty bound- 
ary. Consider the second order hyperbolic initial-boundary value problem 



(<9 t 2 + P(x, D))u = in (0, T) x M, 

u(0, x) = <9 t «(0, x) = for x e M, 

u(t, x) = f(t, x) on (0, T) x dM, 



where 



P{x,D) 



v'det g 
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is a first-order perturbation of the Laplace operator — A 9 on (M,g). Here 
b and q are a covector field and a potential, respectively, in M. We prove 
Holder type stability estimates near generic simple Riemannian metrics 
for the inverse problem of recovering g, b, and q from the Dirichlet-to- 
Neumann(DN) map associated, / — > d v u — i [v, b) g u\gMx[o,T] modulo a 
class of transformations that fixed the DN map. 

1 Introduction and main results 

Let [M , g) be an oriented and compact Riemannian manifold with non-empty 
boundary dM. In this paper we consider the stability of the inverse problem of 
determining a Riemannian metric together with the lower order coefficients of 
the second order hyperbolic initial-boundary value problem [2j from the infor- 
mation that is encoded in the hyperbolic Dirichlet-to- Neumann(DN) map A 9 , 
see ([3| . The question that we want to address is the following: if two hyperbolic 
DN maps are close in an appropriate topology, how close are the Riemannian 
metrics and the lower order coefficients? We now describe in more detailed the 
problem and the results. 

We denote by A g the Laplace-Beltrami. In local coordinates, g(x) = (gij(x)) 
and 




where (g ir ) = (c^) 
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Any second order uniformly elliptic operator with real principal part can be 
written in local coordinates as 



where b is a complex-valued covector field on M and q the complex-valued 
function on M. Moreover, it is self-adjoint w.r.t. L 2 := L 2 (M,dV g ) if and only 
if b and q are real valued, see next section. 

For such P(x, D) let us consider the problem 

( (df + P(x, D))u = in (0, T) x M, 

(2) I u(0,x) = d t u(0,x) = for x € M, 

{ u(t,x) = f{t,x) on (0,T)xdM, 

where / £ Cq(M. + x dM). Denote by v = v(x) the outer unit conormal to 
dM at s € dM, normalized so that g^viVj = 1. The hyperbolic DN map 
A P : H&([0, T] x dM) -> L 2 {[0, T] x dM) is defined by 



Flu 

(3) A P (f):=^-i(u,b) g u 
av 



n 

v l —-\v l b t u 

(0,T)xdM i= i OX 



(0,T)xdM 



where v l = 9 lJ v j ■ This modification of the traditional DN map comes natu- 
rally from the Green identity 

(Pu,v) L 2 {M) - {u,Pv) L *(M) = {^pu,v) L 2( dM) - {u,A P v) L 2 (dM) 

of ([!]). This DN map is invariant under the group of transformations 

(4) g i — ^ :— if* g, b i— > 6* := (p*b — d9 and q t— > q* := (p*q, 

where ip : M — > M is a diffeomorphism with <p\dM =Id and 9 £ C°°(M , K) with 
0| aM = 0, see (19). 

The inverse problem is therefore formulated in the following way: knowing 
Ap, 

1. Can one determine the metric g, the covector field b and the potential q 
up to a transformation that acts on this coefficients by Q? 

2. Could we do this recovery in an stable way? 

The physical interpretation of this problem is to find the speed of the wave 
propagation (i.e, g) and other physical properties (encoded in b and q) inside 
an unknown body by making measurements on the boundary. The initial zero 
conditions in ^ mean that the body is "quiet" initially and we use a disturbance 
/ to recover information about the interior of body; this comes naturally in many 
applications (e.g. geophysics). 

When b — 0, the first question was solved in the Euclidean case (g = e) 
by Rakesh and Symes [TH] using the result of Sylvester and Uhlmann [2H]- For 
the non Euclidean case, the spectral analogue of the inverse boundary value 
problem, was solved for the case b = q = (see also [7]) by Belishev and 
Kurylev [S] using the boundary control method introduced by Belishev [3], [D- 



2 



1 INTRODUCTION AND MAIN RESULTS 



This proof uses in a very essential way a unique continuation principle proven 
by Tataru [29 . Because of the latter, it is unlikely that this method would 
prove Holder type of stability estimates even under geometric and topological 
restrictions. 

For the general case, the spectral analogue of this problem as well as some 
modifications were considered in [T7J, pQ, [TU], [5], [22], Isakov [II], Romanov 
[20] . [15] and references in there. The first question has a positive answer in the 
general case due to Kurylev and Lassas [H] , they used a geometrical formulation 
(see [15]) of the boundary control method introduced in [3]. 

The problem of stability has presented slower progress. Conditional type of 
estimates are typical for such kind of inverse problems. We assume an apriori 
compactness type of condition of boundedness of the C k norm of the metric, the 
covector field and the potential for some large k. This guarantee continuity of 
the inverse map without control of the modulus of continuity. Weak geometric 
conditions of this type were studied in [2J. 

For the Euclidean case and when b = 0, Sun [27] proved continuity of the 
inverse problem. Also for the Euclidean case, Isakov and Sun [T2] proved loga- 
rithmic stability for 2 dimensional domains and Holder stability in 3 dimensions. 
For the anisotropic case b = q = 0, in 23J Stefanov and Uhlmann proved Holder 
stability for metrics close to the euclidean in a C k norm, later they extended 
their result of Holder stability near generic simple metris in [2B] . For the case 
6 = 0, Dos Santos and Bellassoued [5] proved Holder stability for the potential q 
when go is a fixed simple metric, they also recover in a stable way the conormal 
factor of the metric go within the conormal class when q = 0. 

In this paper we prove Holder type of conditional stability when recovering: 
the metric g, the covector field b and the potential q; all at the same time, under 
the hypothesis that the metrics are near a generic simple metric. A Riemannian 
manifold (M, g) is simple if dM is strictly convex and any two point in M can 
be connected by a a single minimizing geodesic depending smoothly on them, 
see Definition Q] for more details. 

Since simple manifolds are diffeomorphic to the unit ball in the Euclidean 
space, from now on, without loss of generality, we consider the case that M = 
0, where f2 is diffeomorphic to a ball in the Euclidean space and has smooth 
boundary. There exist a dense open subset G k {Q), of simple metrics in C k (fl) 
for k ^S> 1 that consists of those metrics for which the x-ray transform is s- 
injective and stable, see definition [3] This set Q k (fl) contains all real analytic 
simple metrics in M, see [25] , 

The main result of this paper reads as follows. Let P be the operator related 
to the metric g, covector field b and potential q as in (jlj, similarly let P be the 
operator related to g, b and q. Suppose that we consider the initial boundary 
problem |2]) for both operators. 

Theorem 1. There exist k 3> 1, < /J, < 1, such that for any go g G k (Q) there 
exist eo such that if 

\\g- 9o\\c(ti) + \\9-9o\\c(a) + \\ b - &llc(fi) < e o 
(5) and 

g,g,b,b,q,q belongs to a bounded set in C k (£l), 
then the following holds. ForT > diam gQ (Q) there exist aC 3 (f2) diffeomorphism 
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ip and smooth function 9 as in Q), such that 

^ \\9 - <?*llc 2 (o) + II& _ ^*llc 1 (n) + Ik - <?*llc(n) 

- C\\h-P - ^-pllfl- 1 ([0,e]xen)^L 2 ([0 : T]xaf2) 

for any < e < T . 

Remark 1. All C fc (f2) norms are related to an arbitrary but fixed choice of 
coordinates. 

The idea of the proof is to divide the recovery in three parts: First we 
prove stability at the boundary using asymptotic solutions pointing in different 
directions to divide the information. Second, we recover the boundary distance 
function from the data and use boundary rigidity estimates in |25j to obtain 
stability for the metric. Third we translate the problem of recovering b and q to 
an x-ray type of problem. We can explicitly recover the x-ray transform along 
geodesies of b and q from the data. We then use estimates obtained in (24] to 
get the desired stability. In each step we must separate the information to get 
estimates for g, b and q separately. 

Acknowledgements: I would like to thanks Plamen Stefanov for suggesting 
the problem and to An Fu for providing his notes on boundary determination. 



2 Preliminaries 

Let us first state some mapping properties for the initial boundary value problem 
([2]) and some remarks about self adjointness. Any second order uniformly elliptic 
operator with real principal part can be written as 

(7) P = -Ag + B + Q 

where B is a complex-valued vector field which in local coordinates has the 
form B = B 3 d/dx 3 and Q is a complex- valued function on M. An straight 
calculation shows that P is self adjoint with respect to L 2 (M, dV g ) if and only 
if 

(8) Re B = and Div B = i 2 Im Q. 
If we write Q as in then we get that 

(9) B = i2fe s and Q = q+ \b\ g + i DivbK 

where (b^y = g^bi. We then see that by |8| the operator |l]) is self adjoint 
w.r.t L 1 if and only if b and q are real valued. 

The direct problem ^ has the following mapping properties, see [IB], [IT] 
and [13]. 

Lemma 2. Let 

(10) F G i x ([0, T];i 2 (M)) and f G Hq([0,T] x dM). 
Then there is a unique solution u{t, x) of the problem 

( {df + P(x,D))u = F{t,x) in (0,T)xM, 

(11) I u(0,x) = d t u(0,x) = for x G M, 
{ u(t,x) = f(t,x) on (0,T)xdM, 
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such that 

u{t,x) G CQO^jff^M)) n C 1 ([0,T];L 2 (M)) 

and 

max {|| M (t)|| ffl(M) + || Ut (t)|| i2(M) } < c(T){ / ||F(t, •)IU=(M)*+||/||^([o,T]xaM)}- 



0<t<T 

'2 







Moreover Af\[ ,T]xdM G ^ 2 ([0,T] x 9M) and 

||A/|[o,T]x3mIU 2 ([0,T]x9M) < c(T){ / -)IIl 2 (M)^ + ||/|lffi([0,T]xaM)}- 

Jo 

To simplify notation we denote from now on, 

(12) II ' II* = II ' \\H^([0,e]xdQ)->-L 2 ([0,T]xdn)- 

where e and T are as in Theorem [T] More precisely, if A is the DN map defined 
above then || ■ ||* is defined as the supremum of ||A||#iqo,t]x<K3) over all / G 
Hq([Q, e] x dfl) such that ||/||Hj([o,e]xan) = 1- This is well defined since one can 
extend / as zero for t > e and this extension of / will be in Hq([0,T] x dfl) 



with f\t=o = 0. Moreover, because of uniqueness for (11), for any < e < T 
and < T < T we have 

(13) ||Ap|| 

Hi ( [0,e] x 9fi)->L 2 ( [0,T'] x dQ.) Hl([0,T] x 9Q)-i-L 2 ([0,T] X 9f2) • 

2.1 Group of transformations that keep invariant the DN- 
map 

We will consider the type of transformations that do not change the DN map. 
Let 

ip : M -> M 

be a diffeomorphism with p\dM =Id, let us denote x = tp(y) then for any 
veC^(M) 

(P{x,D)u,v) L 2 = (<p*P(y,D)(p*u,<p*v) L 2 
where the operatos ip*P is as in with 

(14) p* g , tp* b and ip* q instead of g, b and q 

respectively. Here p* denotes the pull-back with respect to the metric g. Since 
v G C >o (M) is arbitrary we get 

(15) P{x,D)u = Q ip*P(x,D)(p*u = 0. 

Making a change of variable in the DN map operator we see that for any 
/ G C 1 (1R + x dM) 



Ap(/) = Sfe -\v%u | (0T)xaM = -iv%b^p*u 
= (p*v) kd -0-^v) k {p*b) k p*u | v = VKf) 



(0,T)xdM 
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the last equality follows from (15), hence A v ,p = Ap. 

There is also another type of transformation that keep the DN map invariant, 

let 

(16) 9€C°°(M;M),e\ 9M =0. 

We denote Pg the operator as in ([I]) where b is replaced by 

(17) b e = b-d6 

we claim that Ap e = Ap. Similarly as before we can show that 

(18) Pu = P e v = 

where u = e l6 v, then we get that for any / € Cg(K + x dM) since e l6 \gM = 1 



Ap(/) = ^-ii/'M | (QT)x0M = u^-iSb* 
= iu% - §t) v \^ T)xdM = A Pe (f) 



(0,T)xOM 



were the last equation follows by ( 18 ). It is worth mention that transformations 
(14| and (17) commute since (p*b — ip*d9 — ip*b — dip*9. Previous discussion 
justifies the definition of the set of operators that fix Ap as: 

(19) [P(g,b, q )] := {(<p* 9 ,<p*b-d9,<p* q ) : ^S^nd^t"^} 



We emphasis that the recovery of g, b and q is done up to this class (19). 

2.2 Simple metrics and geodesic x-ray transform 

In this section we define the set of generic simple metrics that we will use for 
the proof of Theorem [T] 

Definition 1. We say that g is simple in M, if dM is strictly convex w.r.t 
g, and for each x 6 M, the exponential map exp^, : exp~ 1 (M) — > M is a 
diff comorphism . 

Remark 2. Note that since all requirements for simplicity are open, then a small 
C k (£l) perturbation of a simple metric in O is also simple, so we can extend g 
in a strictly convex neighborhood Hi D il as a simple metric in f2i. 

Consider the Hamiltonian H g (x,£) = (g 1 ^ /2. We denote by (x(t),£(i)) 
the corresponding integral curves of H g on the energy level H g = 1/2. We use 
the following parametrization of those bicharacteristics. Denote 

(20) r_(n x ) :={(z,u)€T*Sl 1 ;z€dSl 1 ,\w\ g = l,(w i u(*))g<Oh 
where v(z) is the outer unit conormal to dQ. Introduce the measure 

dfx(z, u>) — \{oj, v(z)) g \dS z dSu on T_, 

where dS z and dS u are the surface measures on dfl and {w € T*Qi; \w\ g = 1} 
in the metric, respectively. Define (x(t),£(t)) = (x(t; z, cj), £(t; z, lj)) to be the 
bicharacteristics issued from (z,uS) G r_(f2i). 
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For any covector field b = bidx 1 we define the geodesic x-r&y transform I g b 

by 

(21) (Igb)(z,uj) = ( bi(x(t))C(t)dt, {z,u) G r_, 



similarly for any symmetric 2-tensor / = fijdx l dx 3 the geodesic x-ray transform 
I g f, which is a linearization of the boundary rigidity problem, is defined as 

(22) (!,/)(*, w)= J f tJ (x(t))e(tW(t)dt, (z,u)eT_, 

where (x(t),£(t)) — (x(t] z,u),£(t; z,u))) as above is the maximal bicharacteris- 
tics in S7i issued from (z,u) and x % = gij£j = C- 

For a vector field v let d s be the symmetric differential defined by (d s v)ij — 
|(V ' iVj + V,Wj) where Vj are the covariant derivatives. It is known that 

I g d s v = for any vector field v with v\qq = 0. 

Definition 2. We say that I g is s-inyective in f2, if Igf = and / G L 2 (£l) 
imply / = d s v for some vector field v G i?Q(ri). 

Definition 3. Given fc > 2, define ^ fc := G k (n) as the set of all simple C k (n) 
metrics in for with the map I g is s-injective. 

By [25 , for k 3> 1, Q k is an open and dense subset of all simple C k (ft) metrics 
in fi, in particular, all real analytic simple metrics belong to Q k . Also all metrics 
with a small enough bound on the curvature (in particular all negatively curved 
metrics) belong to G k , see [21] and references in there. 

2.3 Interpolation Estimates 

This interpolation result is needed to use the apriori conditions in |5]) . We fix a 
simple metric go G C k for k 3> 1. By remark [2] extend go as a simple metric in 
some ill 3 Let g,g be as in Theorem [l] then there exist A > and £0 -C 1 
such that 

( 23 ) llsllc^n) + Hffllc fc (") - A ' \\9 ~ 9o\\ C (tl) + \\9 ~ 9o\\c(Ti) < £ o- 

The first condition above is a typical compactness condition. Using the inter- 
polation estimate in |30j . 

( 24 ) ll/llaH^^CII/H^II/llc^o) 
where < 8 < 1, t\ > 0, t2 > 0, one gets that 

(25) llff-ffo|| ct( n) <C M 4 fc " t)/fc 

for each £ > 0, if k > t; the same is true for g. For our purposes, it is enough to 
apply ( 24 ) with t, t\ and t% integers only, then ( 24 ) extends to compact manifolds 
with or without boundary. 
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2.4 Boundary normal coordinates and vector field trans- 
formation near the boundary 

In this section we will construct a diffeomorphism <p that fixes the boundary, 
and the smooth function 9 that modifies the covector field, so that near the 
boundary they have convenient properties for calculations. To fix notation we 
sate the following well known proposition about boundary normal coordinates. 

Proposition 3. Let M be a Riemannian manifold with compact boundary dM . 
Then there exist T > and a neighborhood N C M of the boundary dM together 
with a diffeomorphism if : dM x [0, T) — ¥ N such that: (1) <p p (0) = p for all 
p £ dM and (2) the unique unit-speed geodesic normal to dM starting at any 
p € dM is given by 1 1— > f p (t). 

Remark 3. For t > 0, ip p (t) — x so that the distance from x to dM is equal to 
t, i.e. (p p (d(x,dM)) — x for all x in the interior of M. 

For the two metrics g and <?, there exist ipi and <~pi diffeomorphisms like in 
Proposition [3j Then 

(26) p := tpi o ip^ 1 

is a diffeomorphism near dfl fixing <9f2, and mapping the unit speed geodesies 
for g normal to dfl into unit speed geodesies for g normal to <9f2. By [TS] this 
diffeomorphism can be extended to a global diffeomorphism, let us called (p 
again for its extension. 




Figure 1: Boundary normal coordinates 

Notice that g and ip*g have common normal geodesies to <9f2, close to dCl, 
and moreover, if (x' , x n ) are boundary normal coordinates near a fixed boundary 
point for one of those metrics, they are also boundary normal coordinates for 
the other metric, see Figure [T] 

Now let b, b be two covector fields in fl as in Theorem ([!]). Use boundary 
normal coordinates (x',x n ) near the boundary dQ and let b = bidx 1 and b = 
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bidx 1 . Define 

(27) 9(x',x n ) = [ (<p*b n -b n )(x',y)dy. 

Jo 

Extend 9 so that 9 £ C°°(f2; C), notice 9\on = 0. If we use this 9 as in Q, we 
get that near the boundary d£l 

(28) (&.)„ = (ip*b - d6) n = b n . 

Now by the invariance of the DN map under this type of transformations from 
now on we will modify the initial coeficients 

(29) g i-> ip*g, b ^ ip*b - d9 and q ip*q. 

Remark 4. Notice also that by construction of ip = Id+O(eo) m C k ~ 2 ; therefore, 
the metric ip^g also satisfies ([5| with k replaced by k + 3, and for some e' > 0, 
such that £o -> as £q -> 0. The same follow for <£*b — d9 and (fi*q, notice 
that here we are also using that b is close to b as in (I5j). Hence without loss of 
generality we can denote (p*g, (p*b — A9 and (p*q again by g, b and q. From now 
on we will follow this notation. 



3 Stability at the boundary 

We will prove first stability at the boundary following |26j . We consider a 
highly oscillating solution of ([2| asymptotically supported near a single geodesic 
transversal to <9fi. We only need to work locally near a fixed point £o G dfl. 
Let (x',x n ) be the boundary normal coordinates near xq. Consider parameters 
A e K, A > 1 and oj' € K n_1 . Let e and T such that < e < V < T. Fix 
t such that < t < e, and let \ £ C^°(K + x dil) be supported in a small 
enough neighborhood of (i , a;o) of radius e' < min(t , e — t ) and equals to 1 in 
a smaller neighborhood of this point. We define u to be the solution of ^ with 
T instead of T and 

(30) f = e i W-*'-»') X (t,x'). 

One can get an asymtotic expansion for u in a neighborhood [0, T'] x U C K + x O 
of (to, xq) by looking for u of the form 

N 

(31) u = e iX ( t -* {x ' ul))y £ j \- : >A j {t,x : ui) + 0(\- N - 1 ) 

3=0 

in CH^T'JjL^t/)), where TV > is fixed and u = (uj',uj n ) is such that 
^ g lJ (xo)uJiUJj — 1, ^2 9 lJ (xo)vi(xo)ujj < 0. In U the phase function solves the 
eikonal equation 

with the extra condition < it is uniquely solvable near xq. In our 

coordinates, the metric g satisfies gi n = g™ = 5i n for i — 1, . . . , n and d/dv = 
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—d/dx n . Notice that (32) determines 

(33) io n (x) = > 0,Va; G Undft. 

In [0, T'} x U the principal part Ao of the amplitude solves the first transport 
equation. 

(34) LA o = 0, A Q | K „ =0 = X 
and the lower order terms solve the transport equation 

(35) iLAj = -(<9 t 2 + P)Aj-x, A i U„ =0 = 0, j>l, 
where 

(36, t ^ l+1 E.'Si-»EA^ v . 

The construction of u guarantees that A,-, j = 1, . . . , N are supported in a 
small neighborhood, depending on the size of sup x> of the characteristics issued 
from (to,Xo) in the codirection (l,w). Therefore by the way we choose e', the 
term 

JV 
3=0 

in (31 1 satisfies the zero initial condition in Moreover, ujy satisfies the 

boundary condition Ujsr = f with / as in ( |30| , provided that T' is such that 
< T' — t , is small enough. Write 

U = UN + W. 

Then w = w t = for f = 0, and tu|(o,T')xdn = and (df + P)w = 
e iA(t-0) A -w^2 + p}A N . By Lemma we obtain ^ with the estimate 
remainder in the C 1 ([0, T']; L 2 (U)) norm. We emphasize that it is important 
that T' — to is small enough so that the wave does not meet dVl again (if it does, 
we need to reflect it off the boundary, as in next section). 

In order to get stability at the boundary we are going to take advantage 
of the freedom that we have of choosing the initial codirection cj G R™ of the 
solution u. For that we need the following lemma, we postpone its proof to the 
end of the section. 

Lemma 4. Let M be a Riemannian n-compact manifold with continuous metric 
g. For any xq G M and any chart (W, ip) containing xq there exist U C W open 
neighborhood of xq and u>i(x), . . . ,u)n(x) G S^M = {£ G T X M : |£| s = r} for all 
x G U, such that any symmetric 2-tensor field h(x), and any 1 -tens or field j3(x) 
there exist a constant C > such that 

max \hij(x)\ + \f3j(x)\ < C max \T uli [h(x), /3(x)]\ for all x G V 

l<i,j<n l<i<N 

where N — n(n + 2)/2 and T u t x )[h(x), /3(x)] — hij{x)uj l {x)uj 3 (x) + ^(x)^ (x). 
The main result of this section is the following: 
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Theorem 5. For any fj, < 1, m > 0, there exist k 3> 1, such that for any A > 0, 
if 

g,g,b,b,q,q are bounded in C7 fe (f2) by A, 
then there 3C > depending in A and Q, such that 

ft) su Pxe90 ,| 7 |< m \8H{g - <fg)\ < C\\A P - Ap^ 2 ™ 

(ii) sup xGan , |7| < m \ffr(b - &*b - dff))\ < C\\A P - Ap||^ /2m+1 



(Hi) sup^gan. 



h\<r 



\drr(q-<p*q)\<C\\A P -A p \ 



where ip and 8 are as in (26) and (21), respectively. C k norms are taken as in 



Remark^ 

Proof. By Remark^ without loss of generality, we denote f*g, (p*b — d9 and 
ip*q again by g, b and q. To simplify notation let S — ||Ap — Ap||*. In this proof 
we will denote C by various constants depending only on f2, A and the choice 
of k S> 1 in Theorem [T] Using local boundary normal coordinates as before, let 
I x V C M+ x dfl be a neighborhood of t = t , (x',x n ) — (0,0), where X = 1. 
Since = — gfrr and by Lemma [SJ we have 



i\(t—x-uj) 



(37) A P f = e 



in L 2 (I x V); similarly for Apf. Then 




- - i&„) Aj j + OiX-"- 1 ) 



d<j> d(j) 



1 



dx n dx n i\ 



,— i\(t— x-ui) 



(A P f -A p f 



8A l _dA l 



dx r ' 



in L?(l x V). Hence, 



d(> d<p 
dx n dx n 



dx n 



C 



i\ 'lb, A, /,,,!,: +0(\- N - 2 ) 



C 



L 2 (V) 



< V (^l|/||HH[0,e]x9n); . ^ 



Notice that ||/||H 1 ([o,£]x9n) < CA. Now take A — > oo above to get 



(38) 



dx n dx n 



< C8. 



By the eikonal equation ((32]) , in V C 9fi, we have 
(39) 



1/2 



a, (3=1 



and similarly for ui n = d<p/dx n . Since ||(w ra ) 2 — (oj n ) 2 \\L 2 (V) = 0(£) by (38) 
Then we get 

\\{g<*P{x)-~g^{%))uc.MW { v)=0{8). 



11 



3 STABILITY AT THE BOUNDARY 



and by Lemma [4] we have 

(40) \\g-~g\\v(v)<C5. 



Using compactness the manifold we extend (40) to the whole dfl hence: 
(41) \\g - g\\man) < CS. 

We use interpolation estimates in Sobolev spaces and Sobolev embedding theo- 
rems, to get that for any m > and \i < 1 



(42) 



\\g-g\\c-(on)<CS^ 



provided that k 3> 1. 

To estimate the difference of the first normal derivatives o f q a nd g and the 
difference between b and b we use up to the principal part in ( 37 1 , as before 

- g - , M - M„) - -.-•«" ->(M«) - a,(«» + a (£ - g) 



3=1 



j = l 



to obtain 
(43) 



dA dA 



^-i(b n A -~b n A ) ^CiXS + S + X- 1 ) 



The r.h.s above is minimized when A = 8 1 ^ 2 
dA Q dA 



(44) 



i(b n A - b n A ) 



L 2 (V) 



dx n dx 

The transport equation on (to — ei,to + £i) x V implies 

dA 



< cs 1 / 2 . 



2w, 



dx r ' 



+ A g( f>~2iY^ !l' JI >.~, = 0. 

i,j=0 



Assuming the w„ , uj n > r/ > for small r/ then by (32) and since g™ = Si n , 
o o- f, w « 9det. 9 <9 2 </> a 



w n 9 det g 1 dg 



a/3 



2 det g dx 11 2oj n dx 7 



-uj a uj/3 — 2ilfwp + R, 



where R involves tangential derivatives of g that we can estimate by ( 41 ) and 
d(f)/dx n that we can estimate by (38). Therefore by (44), 

u! n ddetg uj n ddetg 
2 det g dx 11 2detg dx n 

(45) 



- + 



1 dg a >P 1 dg a >P 

~2i{b Lo p - 6%) = OiS 1 / 2 ) 
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in L 2 (V) for all oj's as above. Setting <J = 0, we get 



1 d dct g Id det j 



2&etg dx n 2 det g dx r 



= 0(«5 1/2 ), 



and then since we can estimate the difference of the metric by (41) we obtain 
d det g d det g 



dx r ' 



dx r ' 



0(5 



l/2\ 



LHV) 



This last equation together with ( 45 ) implies 

1 dg a # 1 dg 01 ' 



2u n dx r ' 



J a 0J/3 



2uj n dx r 



0{8 



1/2-1 



L 2 (V) 



since ui n , u) n > r\ > and using again ( 38 ) we have 

~dx~ n dx^~ J WqW/3 ~ ^ n ~ 4u}nlb 

Let us now take uj 2 = = 1/2, then 



0(* 1/a ), 



L=(V) 



and also w' belongs to {to' € T x >dfl : \u)'\ g i — 1} where 5' is the induced metric 
of either g or g to 50. Using now Lemma [4] and the fact that near the boundary 
b n — b n = 0, see (28), we have 





\\b-b\\mv) 


+ 


dg 

dx n 


dg 

dx n 


< cs 1 / 2 . 

L 2 (V) 


Again by compactness we get, 








(46) 


\b- £>||z2(0fi) 


+ 


dg 

dx n 


dg 

dx n 


< C5 1 ' 2 . 

L 2 (dQ) 


As before using interpolation 


and Sobolev embeddings theorems 


(47) 






+ 

(9J1) 


b-b 


c m (dn) 



for any m > and fi < 1 as long as k 3> 1. 

To estimate the difference of the second normal derivatives of q and g, first 
normal derivatives of b and b and the difference of q and q we use ( 37 ) up to the 
A -1 to get 



dx n dx n 
Choose A = <5 -1 / 4 to obtain 

(48) 



L 2 (V) 





dA l 


dx n 


dx n 
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3 STABILITY AT THE BOUNDARY 



Using the equation iLAj = — (df + P)Aq restricted to I x V we get 
dAi 1 ddetgdAo d 2 A 



2iuj rl 



dx 11 2 det g dx n dx n d 2 x n 



q, 



this together with (|38|), (|4l|) and (|48|) gives 
(49) 



d 2 A n _ d 2 A n „„ . 
2w n ^- - - 2w„^— - + 2w„g - 2uj n q 



Q2 x n 



d 2 x n 



= 0{5 1 ' 4 ). 



L 2 (V) 



Assuming u n ,Q n > r\ > for small r\ and taking normal derivatives in (34) and 
restricting it to / x V we have 



2u> n 



d 2 A 

d 2 x n 



u n d 2 det g d 3 ^ a adb a 

219 d^^ + R 



+ 



2 det g d 2 x n d 3 x n 

u n cfdetg 1 3 2 g^ af3 db a 
2detg Wx" + 2^9^ Wq ^~ 215 dx-^ +R 



where R consist in term that we can estimate by ( 41 1 and ( 46 ) . Hence we get 
uj n d 2 det g Q n d 2 det g 



2detg d 2 x n 2&etg d 2 x r 



+ 2u n q - 2cu n q+ 

1 d 2 g a ^ 
2uj n o^x™ 2oj n a z x n 



1 d 2 g a >P 
, db n 



in L 2 (V). Again, setting u/ = we have 
(50) 

Now since w n — w n = 0(8) then 



1 d 2 det g Id 2 det g 
2detg d 2 x n 2det<? d 2 x n 



2q-2q 



L 2 (V) 



Oid 1 / 4 ). 



1 

2uj,, 



Q2 a ,f) Q2~ g »,P\ f dh 



d 2 x T 



d 2 x n 



dx r 



5*0 9b a 

y dx n 



LU/3 



0{5 



1/4N 



L 2 (V) 



We use similar reasoning as before. First we use Lemma [4] and then compact 



ness to get that 



db _ db 

dx 71 dx n 



together with ( 50 ) gives 

(51) \\q- q\\v>(en) + 
As before we get 

(52) \\q- q\\c^(dfi), 



L 2 (dU) 



and 



8? g _ _&g_ 

d 2 x n d 2 x n 



db db 






+ 


dx n dx n 




L 2 (dQ) 



L 2 (dfl) 

d 2 g d 2 g 



are 0(J 1/4 ), this 



db 


db 




dx 11 


dx n 


1 









d 2 x n d 2 x n 



d 2 g d 2 g 



L 2 (dU) 



d 2 x n d 2 x n 



< C*<W 4 . 



for m > and ji < 1. Proceeding by induction we prove the theorem. 



□ 
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Proof of Lemma [JJ First notice that by rescaling it is enough to proof the the- 
orem for r = 1. Let ^i(xq), ■ ■ ■ , £iv( x o) be like in Lemma 3.3 in [5] related to 
x . Consider 

= jH X °l Vi = 1, . . . ,n and m = 1, . . . , N. 

By continuity of the metric, for any e > 0, there exist a small neighborhood U of 
x were |£ m (xo)| a > 1/2 and sup lGy max ra , i |^(a;) - C m ( x o)\ < £• Now notice 
that in Lemma 3.3 in [8] the determination of [ft, 0\ is done by inverting any of 
the linear transformation Tq : [ft, j3] H> (T Wl [ft, /3], . . . , T WJV [ft,, /?]), whose inverse 
is also linear. Choose e small enough so that linear independence is preserved, 
then we can take C — sup xg(7 1 | ■ were U might be a smaller neighborhood 
of x . □ 



4 Interior Stability 

We first estimate the difference in the metrics and its derivatives following the 
argument in that is based on the boundary rigidity result in[5S]. We then es- 
timate the difference in the co- vector fields and the potentials using semigeodesi- 
cal coordinates related to a point zq (£ but close enough so that simplicity 
assumption is still valid. We use the phase function 

4> — X n — d(x, Zq). 



Then the transport operator ( 36 ) becomes 

„ d „ d „ , 1 9det q 

(53) L - 2 m +2 o^~ 2 ^ + vm^ 

and can be solved explicitly. We get non-weighted integral transforms of b and 
q. We then apply Holder stability theorem in [21] for vector fields and functions 
to prove the theorem. 

For technical reasons we will need to make a 5-small perturbation of g, b 
and q so that near the boundary they coincide with g, b and q. We will follow 
notation as in Remark [4] 

By Theorem[5j one has that for any m > 0, there exist < /i < 1 and k 3> 1, 
such that 

(54) sup \IF(g-~g)\<C5». 

Let m > be any integer and let \ € C°°(M), x = 1 for t < 1 and \ — for 
t > 2. Set 

(55) ~g 1 ^g + X {s- 1 / AI p(x,dn))(g-~g), 

where M = 2m /n Using a finite Taylor expansion of g — g up to 0((x n ) M ), 
where x n = p(x,dil) and estimate (54 1, we see that 

(56) \\~g 1 -~g\\cm(n ) <C8 fl - m/M = C6»/ 2 
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Is important to notice that g\ depends on m, but we will only need (55) for 



large fixed m. In particular the estimate above implies that g\ is also simple for 



S -C 1. As before without loss of generality we can assume that (23) are still 
true for g and g\. We extend g and §\ in the same way as simple metrics in a 
neighborhood Cli D CI. The advantage we have now is that 

(57) g = g~i for - 1/C < p(x, dCl) < 6 1/M . 

and hence by strictly convexity of the boundary there exist a constant C K such 
that 

(58) ( Pg -paJlaaxti = if p g (x,y) < C K 5 X ^ 2M \ 



where p g denotes the distance function with respect to g. Moreover, using (56) 
we obtain 

(59) \p s (x, y) - Pgi (x, y)\ < CS^ 2 Vz, yen. 

We define similarly 

bi = b + X (S- 1/M p(x, dCl)) (b - b) 

(so) ; ; 

q 1 =q + X (S- 1/M p(x,dny) (q-q) 
and use Theorem [5] to get that 

(61) Wh - b\\ cmm + Wq, - q\\ cmm < CS 1 ^ 

and 
(62) 

We proceed to the proof of the Main Theorem: 



b = b 1 &ndq = q-y for - 1/C < p(x, dfl) < S 1 ^ 1 



Proof of Theorem^ Recall the notation S = ||A — A||*. It is enough to prove 
the proposition for 5 <C 1. We use notation as in Remark |4j In what follows 
we denote p < 1 constants arbitrarily close to 1 that might change from step to 
step. We also denote by C various constants depending only in ft, A and the 
choice of k in Theorem [TJ 



We first modify the mertic g, b and q by g\ , b\ and q~i as in ( 55 1 and ( 60 ) . We 
denote this extensions again by the same notation. From now on objects below 
related to g are without tildes and those related to <?i are with a tilde above (no 
subscript 1). Fix xo,yo £ dfl. We will construct an oscillating solution related 
to g similarly to the one used in the previous section, but in this case we want it 
to go all the way from xq to yo. Consider the geodesic from xq to yo, extended 
from CI to z £ Cli such that the geodesic segment [zoi^o] C fii\r2, see Figure 
[2j Assume also that p g (zo,x ) > 1/C > 0. 

Consider (x' , x n ) global semi-geodesic coordinates related to zo, where x n is 
the distance from zq to the point p and x' is a parametrization of the angular 
variable so that x' = const, are geodesies issued from z$ with x n as the arc- 
lenght parameter, as in Lemma 4.2 in [25] . We have in this coordinates 

x n = p g (p, z ) and g in = S in . 
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Figure 2: Solutions near semigeodesical coordinates. 



The prove is divided in three step: 
(a) Metric Stability: We follow and use boundary rigidity estimates in 
Some simplifications to the argument are presented. We claim that 



(63) 



I \Pg - Pgi I lc(anxdQ) < OS* for < fj, < 1. 



Assuming the claim and using estimate ( 59 ) we have 

\\p g - PgWc(dQxdQ) < for < \i < 1. 
We then apply estimate in Theorem 1.8 in [25] to get 
(64) \\g-g\\a»m < Ctf M for </*<!. 



We will sketch the proof of the claim ( 63 1 since is similar to that in [5^1 . We write 
the explicit constr uct ion as in th ere since it will be used in part (b) and (c) of the 
proof. In view of (58), to prove (63) we assume w.l.o.g. that p(xo,yo) > CW • 
We claim that 



(65) 



\p 2 Jx ,y Q )-p 2 ,(x Q ,y Q )\ =0(6") 



with some < fi < 1 uniformly w.r.t. xo and yo, which implies (63). 

Let T > diamfi, w.l.o.g. assume that T — diamfi is small, see (13 1. For 
< t < 1 define 



(66) 



U = £l|X dfi; \t - t 1 +p(x,a;o) < 



where Ci 3> 1 will be specified later. Choose a cut-off function < \ < 1?X S 
Cq°(1R+ x 957) such that supx C f, and x = 1 in a set defined as U but with 
C\ replaced by 2C±. One can arrange that \d t x\ + |Vx| < C<5 _1 /( 2M ). 

Set </>(.t) = z ). Then, by simplicity assumption, since z Q € fi^, we have 
that (f> E C 1 (0), and <fi solves que eikonal equation 



9x J to? 



1. 
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Now we construct a solution u of ([2]) for 

/ = e -iA( t -0) x . 

There is a need of reflecting the solution ones it reaches the boundary to get 
the zero boundary condition, so the solution u is the sum of the projected wave 
and the reflected wave, u = up + ur with 

ui(t, x; A) = e - 1A (*-*( x )) (A (t, x) + A _1 Ai(t, x) + fl(t, x; A)), 
u B (i, as; A) = e -**(*-*(*)) (i ( i; x ) + \- 1 A 1 {t 1 x) + R(t, x; A)). 

where 

C 

(67) ll-R||c 1 ([o,T];L2(n)) + H-Rllc^io.r];^^)) < 
This solutions satisfy the following equations: 

(68) LA = 0, A \ u = x; i LA X = ~{d 2 t + P)A , A j \ u = 



(69) LA Q = 0, A Q \ V = -A \v; iLAt = -{d? + P)A Q , A x \ v = -A x \ v 

where V C K+ x 9il is the image of U under translations by all geodesies issued 
from Zq and passing through U. The phase function <f) still solves the eikonal 
equation with boundary condition <j>\y = (f> and is unique solution with gradient 
pointing towards the interior of f2 (the opposite solution in <fi). All this boundary 
conditions are assumed to be extended as zero in the rest of the boundary. 

We can solve the solve the transport equations ( |68| and ( [69] ) in a neighbor- 
hood of the geodesic connecting xq and j/o of size 0(lr' ( 2M >) , and by simplicity 
assumption, this solutions can be extended all the way to V. If C\ in (|66|)is 
large enough, then U and V are disjoint sets 

Because of the strict convexity of dfl, each component is of size Of^S 1 ^ 2 ^), 
at a distance bounded from below by the same quantity by assumption. Denote 
by B(y,r) the ball centered at y with radius r. Then V contains the set Vq = 
ynBfeo.^VC ), such that on V , we have A > l/C > 0. Above, C° 
is chosen so that Vb is contained in the translation of the set {\ — 1} under 
geodesies issued from zq. 

Next, we construct a similar solution u related to g. We construct first a 
phase function <f> as <fi = p{x, zq). It solves the eikonal equation 



(70) 



E 



dx % dxi 



4>\u = <t> 



The later equality follows from (57). The other properties of u are similar to 
those of u. Let Vq be defined as above, but associated to g. 

We claim that Vq PI Vq ^ arguing by contradiction. Notice that 



(71) 



Wo 



, W <^" A2 ^A 2 + ^) 
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where the first constant is independent of S (but it depends on A and eo in 
(23 1). To prove this, we only need to estimate AqW4> on U, and use that 
\dtx\ + |Vx| < CS-WW. On V , we have 



(72) 



A P f = -2iAe iA( ^ f) A 



d<\> 
dx n 



O s (l) 



in L 2 (Vq) and similarly for Apf. Notice that 



(73) 



dx< 



> ,5«/(2M) j C 



L 2 (V ) 



because area (Vb) > S n ^ 2M '/C. Since Vb and Vq do not intersect and by (71) 
we obtain 

X 2 S n/(2M )/c _ C{S) < ||Ap/ _ Ap f\\l 2[0tT]xm < CA2(5 1W(2M) + ^ 

Dividing last inequality by A 2 and taking A — > oo we get a contradiction. H ence 
5 £ Vb n Vb, using simplicity assumption and triangular inequality we get (65) 
and this concludes section (a). 

(b ) Magnetic Field Stability: For this part we will use the stability of the princi- 
pal part in the solution constructed in (a) and stability of the 1-tensor geodesic 
x-ray transform. We will use sharp estimates of the 1-tensor x-ray transform ob- 
tained in |24j . Stability for this 1-tensor geodesic x-ray transform was previously 
know, see [21] and references in there. As before 



\w-([a,T\x.dn) - + C(S) 
idependent of 5, but it depe: 

\\A P f-Apf\\ 2 W]xdn <C\ 2 5 + C(6). 



where the first constant is independent of S, but it depends on A and eo in ( 23 1 
We also have 

|2 ^- si\2 



Using this two last equations, (72) and part (a) we obtain 



\A -A \\ L2{Von y o) <CS + C(S)/X, 



taking A — > oo we get 
(74) 



\A -A \ 



L2(v nv ) 



< CS 



In this coordinate system, remember 
comes 



x n , the transport equation (36) be- 



(75) 



d 



d 



1 



d dct g 



dt dx n 2^/detg dx r ' 



- ib n )A = 



with same initial conditions as in ( 68 1 . Using the method of characteristics or 



a change of variables we can compute explicitly this solution and get 

If we stay near t = x n , then x = 1 so 
(76) A (t lX )=^ 



' detg(a;',0) \ 1/4 i /*" b n (s)ds 
det g{x) 1 
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Using ( |74| ), and ^64J) we get that 

(77) ||gi Jo 1 " bn(s)-6„(0,s)ds 



^Iz^VbnVb) - C6>1 - 



Remember that we modified the covector field to that b— b = in a neighborhood 
of the boundary containing Sli \ f2, then b — b belongs to C 1 (J7) and by ^ and 
(61 ) we have 

\\b-b\\cHn) <C(6 + W 2 ). 
Since i5< 1 we can use a Taylor expansion of e tx near zero to get 



(78) 



b n (0, s) - b n (0, s)ds 



< C5". 



L 2 (v nv ) 



This is the coordinate representation of the x-ray transform algon the geodesic 
starting at z and going all the way to uq. Until know we had a fixed coordinate 
system associated to zq. Since all constants C are uniform with respecto to xq 
and yo We can then shoot rays in all directions and we can move z$ around 
all 90x- We note that since we modify the covector field near the boundary of 
£1 there are no non-zero tangential rays that are being integrated over. Hence, 
since I g : £ 2 (f2i) — > L 2 (T_(n,i);dfj,) is bounded by [H] we have 



(79) 



\\I g (b-b)\\ LHr _ {ni) , d ^ <C6». 



Using interpolation and the compactness assumption of the metric and the po- 
tential we can estimate 



By Theorem 4 in 



\\I g (b-b)\\ CH r-(n M <C5^. 
we know 



(so) \\b-b x \\ L . {n) <c\\i;i g (b-b)\\ m{n) 

with 

were y X £ denotes the maximal geodesic in fii that passes through x with codi- 
rection £. Hence 

(81) \\b-b\\ L 2 (ni) <C6^ 



for some < /i < 1. Now by (61 ) and since b = b\ on Q,\ \ £1 then 

(82) \\b-b\\ LH n)<C6». 
By interpolation and compactness we get 

(83) \\b-b\\ cl(n) <CS^. 

(b) Potential Stability: Finally for the potential we use the next term in the 
expansion of the previous solution and stability estimates for the geodesic x-ray 
transform of functions in [23]. We have that the DN-map is given by 

Ap/ = e- iA(t -« (-2i\^A -2i^A l + {d v A -d u A Q )\+O{\- 1 ) 
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in L 2 (Vo) and similarly for Ap. Notice that since the metrics g and g are equal 
in a neighborhood of the boundary then all rays can be taken transversal to the 
boundary. Hence 

d , . d , . d 



dv 



with \a(x)\ > 1/C in Vo fl Vo , where d v ± is tangential to the boundary. Since 
we can estimate tangential derivatives we get 

-2i^(A 1 -A 1 )+2a(x)^-(A -A ) = Apf-A p f+2ia(x)\(A -A )+O(X- 1 )+O(5^ 
av ox 



in L 2 (Vq D Vo). Now using the explicit form of the amplitude ( 76 1 to estimmate 
A — Aq we obtain 



-i^{B x -B x ] 
ov 



< C(X6 + l/X) + 0(5"), 



L 2 (v nvb) 



taking A = S x l 2 we get 

(84) 

and 



A x - A y 



L 2 (V) 



(85) 



d d 



1 d det g 



ib n A, = (d 2 + P)A . 



^dt dx n 2^/detg dx 11 
By ^ then 

(d 2 + P)A a = d 2 A a - A g A Q - qA - \b\ 2 g A - i (2(6, d4>>» + (Vivb)A ) . 
Again, near t = x n we get 



M*>t) = WT) / [v(x',s)A g A - (detgix*)) 1 ^- (det 5 (x')) 1/4 |6| s 



where 



+ i 2r)(x', s) (2(6, dA Q ) g + (Divb)A )]ds 



t]((x, 0)', x n ) = (det gix)) 1 '^-' 1 /o"" b »C s ) ds 



Now by ( 84 ) and previous estimates we have 



since we have estimates for the metrics and the covector field this implies 



q n {s) - q n (s)ds 



< C6". 



L*{V) 
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Comparing this with (78 1, we ague as before, we obtain an invariant represen- 
tation of the previous inequality 

Wlgiq-fillL^r-^dri < Ca- 
using interpolation and the compactness assumption of the metric and the po- 
tential we get estimate (remember we changed the potencial so that q = q near 
the boundary) 

\\I g (q - g)llci(r_(n l)id „) <C5». 

We then apply an stability estimate for the geodesic x-ray transform as in The- 
orem 3 in [24) and interpolation to get 

H?-g1lc(n)<Ctf" 

for some < \i < 1. □ 
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